Can antiferromagnetism and superconductivity coexist in the high-field paramagnetic 

superconductor Nd(0,F)FeAs? 
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We present measurements of the temperature and field dependencies of the magnetization 
M(T,H) of Nd(Oo.89Fo.n)FeAs at fields up to 33T, which show that superconductivity with the 
critical temperature T c ~ 51K cannot coexist with antiferromagnetic ordering. Although M(T, H) 
at 55 < T < 140K exhibits a clear Curie- Weiss temperature dependence corresponding to the Neel 
temperature Tjv « 11 — 12K, the behavior of M(T, H) below T c is only consistent with either param- 
agnetism of weakly interacting magnetic moments or a spin glass state. We suggest that the anoma- 
lous magnetic behavior of an unusual high-field paramagnetic superconductor Nd(Oi_ a; F :E )FeAs is 
mostly determined by the magnetic Nd ions. 

PACS numbers: 74.20.De, 74. 20. Hi, 74.60.-w 



The recently discovered superconductivity in the fam- 
ily of doped oxypnictides [l[ exhibits an unusual coexis- 
tence of the rather high critical temperature T c ~ 26 — 57 
K with strong antiferromagnetic (AF) correlations and 
apparent magnetic ordering. Indeed, neutron 0,1, pSR 
[j, HI, NMR @, 0| and Mossbauer d, H measurements 
have revealed complex magnetic structures including a 
spin density wave AF instability at Tjv — 140 — 150 K 
accompanied by the resistivity drop and monoclinic dis- 
tortions of the tetragonal lattice of the parent semimetal 
compound (Io| . Doping with F or O reduces Tn and 
seemingly eliminates manifestations of the AF ordering 
in the temperature dependence of the resistivity p(T) in 
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This be- 



oxypnictides with the optimum T c 
havior poses a fundamental question if the AF order- 
ing is merely a consequence of the nesting of the Fermi 
surface in the parent semimetal or it does play a ma- 
jor role in superconductivity. Recent theoretical models 
have proposed that magnetic correlations are instrumen- 
tal for the unconventional superconductivity in oxypnic- 
tides [II [II El, I3i Ell i which may also be important for 
understanding their extremely high upper critical fields 

In this Letter we address the specific form in 
which the magnetic correlations manifest themselves 
both in the normal and superconducting state of the 
Nd(Oi_ 2: F :E )FeAs using high magnetic fields H up to 
33T capable of destroying any global AF state with 
Tn < iiBH/ks — 22 K, where ps is the Bohr magne- 
ton. Our measurements of the magnetization M(T, H) 
in a wide range of fields < H < 33 T and temper- 
atures 5 < T < 140 K have shown a strong param- 
agnetic contribution in M(T, H) both at T > T c and 
T < T c . In the normal state M(T) exhibits a clear 
Curie- Weiss temperature dependence with an apparent 
Neel temperature Tjv — 11 — 12 K, which might sug- 



gest a coexistence of superconductivity and antiferro- 
magnetism. However, our data at T < T c show instead 
that M(T, H) in the superconducting state is inconsistent 
with a long-range AF ordering, indicating that supercon- 
ductivity effectively cuts off the AF interaction between 
magnetic moments evident at T > T c . We suggest that 
Nd(Oi — kFk)FcAs may be regarded as a paramagnetic 
superconductor with unusual magnetic behavior differ- 
ent not only from conventional superconductors but also 
from La(0i_ 2; F a ;)FeAs oxypnictides. This difference may 
result from a large magnetic moment p = 3.62/!^ of a free 
Nd 3+ ion (p = 3.58/is for Pr 3+ ), as compared to p = 
for La 3+ and p = 0.8p B for Sm 3+ Thus, m(T,H) 

of Nd(0,F)FeAs is mostly dominated by the Curie- Weiss 
behavior of Nd 3+ , while contributions from Fe 3+ with 
p = 0.25 - 0.35/xb d H and temperature-independent 
van-Vleck and Pauli paramagnetism become negligible. 

The Nd ( Op . 89 Fp. 1 1 ) Fe As poly crystalline samples were 
made as described in Ref. [l3[ ■ The sample with approx- 
imate dimensions of 1.4 x 1 x 0.5 mm 3 , volume V" w 0.7 
mm and weight s» 4.5 mg had p(T c ) ~ 0.53 mOcm, 
T c w 51 K, and about 90% of the theoretical density, 
7.213 g/cm 3 . The magnetic susceptibility %(T, H) and 
magnetic moment M(T,H)V were measured by SQUID 
magnetometer, Vibrating Sample Magnetometer (VSM) 
in a 14T superconducting magnet and a resistive 33T 
magnet at NHMFL. Shown in Figs. 1 and 2 are the 
zero field cooled (ZFC) and field cooled (FC) M(T) 
curves, which change radically upon increasing field. For 
H = lOmT, M(T) exhibits a standard superconducting 
behavior with the diamagnetic onset close to T c . How- 
ever, increasing field to a rather modest value of 0.5T 
(given H c2 (0) - 50-100Tfor oxypnictides [II])), shrinks 
the irreversible diamagnetic loop, while significantly in- 
creasing the paramagnetic component of M(H). Further 
increase of H to 1.5T, practically eliminates the differ- 
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FIG. 1: ZFC (black full dots) and FC (red open dots) M(T) 
at: 0.1, 5 and 15 kOe. Inset shows FC M{T) at 0.1 kOe. 
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FIG. 2: ZFC magnetization at different fields. The dash line 
shows T m extracted from the data in Fig. 3. 

ence between the ZFC and FC curves, both exhibiting 
a strong paramagnetic behavior. Fig. 2 shows in more 
detail how this field-induced superconducting to param- 
agnetic transition develops. 

To gain a further insight into the nature of the param- 
agnetic signal not masked by superconductivity, M(T, H) 
between 55 and 140 K was measured by a SQUID mag- 
netometer. The results shown in Fig. 3 unambiguously 
indicate the AF Curie- Weiss behavior, M _1 oc T + T m 
with T m 11 — 12K being practically independent of 

at < B < 5 T. This may indicate AF ordering at 
the Neel temperature Tn — T m , below the superconduct- 
ing transition and a coexistence of superconductivity and 
antiferromagnetism in a variety of forms, for example, a 
quantum critical point as a function of doping. To see if 
this is indeed the case, we measured M(T, H) below T c at 
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FIG. 3: The AF Curie- Weiss linear temperature dependence 
of M _1 (T,.H") (x T + T m in the normal state. The linear 
extrapolations defines T m « 11 — 12K. The inset shows the 
collapse of the inverse susceptibilities x _1 (T, H) = H/M for 
different fields < H < 5 T onto a single straight line. 



very high magnetic fields, which enabled us to suppress 
the irreversible magnetization due to pinned vortices and 
to probe M(T, H) in the superconducting state at fields 
H > Tn/h sufficient to destroy the global AF state. 

Shown in Fig. 4a is the M(H) loop at 4.2 K for 
the field perpendicular to the broader sample face. One 
can clearly see a superconducting hysteresis loop on top 
of a strongly paramagnetic, field-dependent background. 
Taking the mean value of the two branches of m(B) as the 
paramagnetic magnetization, M p = (Mj + MjJ/2 plus 
the small reversible magnetization of the vortex lattice, 
we define the magnetization AM = M-f — M± due to the 
critical state of pinned vortices. Using the Bean model, 
we estimate the critical current density J c (4.2if) ~ 10 
kA/cm 2 at self field. M(H) in Fig. 4a increases with 
field and tends to level off at higher H. To probe the re- 
gion of higher fields, M(H, T) was measured in a 33T 
resistive magnet. The results presented in Fig. 4b, show 
that M (4.2K, H) does saturate, while M(T, H) at higher 
temperatures increases much slower as H increases. The 
width of superconducting loop decreases as T and H in- 
crease, consistent with the usual behavior of J C {H, T). 

To interpret our data we used the mean-field AF the- 
ory in which two magnetic sublattices with the magneti- 
zations Mi and Mi are controlled by the molecular fields 
Hi = H + aMi - bM 2 and H 2 = H + aM 2 - bMi, 
where a and b quantify the exchange interaction with 
nearest neighbors of the same spin sublattice (a) and 
the sublattice with the opposite spins (b) [2l[. The self- 
consistency equations, 2M, = M Bj(^Hi/T), i = 1,2 
determine the total magnetization, M(T, H) — Mi + M 2 , 
and the AF order parameter L(T, H) — M\ — M 2 where 
Mq = njj,, ji is an elementary magnetic moment, n is 
their total density, Bj{z) = [(2J + 1)/2J] coth[(2J + 
l)z/2J] — (1/2J) coth(z/2J) is the Brillouin function, and 
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FIG. 4: (a) Low-field M(H) loop measured by VSM at 4.2 K. 
The red line shows M = Mo tanh(/uB/T) with B = H + AtyM, 
M = fin for a = 0.69n B , M V = 0.135 emu and n = 3 x 10 22 
cm -3 . Inset shows m(H) and 1(H) calculated from Eqs. ([2j- 
(J6]) for 7 = 0.5, and m = tanh(/i) (dashed line). Here the left 
panel shows the second order transition for 6 — 5/11, and the 
right panel shows the first order metamagnetic transition for 
= 0.1. (b) High-field M(H) loops at different T. The red 
lines show M = M tanh(^B/T) with u = 0.67/ib, MoF = 
0.133 emu, n = 3.1 X 10 22 cr 
0.123 emu, n = 2.25 x 10 22 cm- 
MoT/ = 0.112 emu, n = 2.05 x 10 22 



~ 3 (4.2K); fi = 0.85ms, M V = 
~ 3 (12.5 K); fj, = 0.85u B , 
cm" 3 (20 K). 



J is the maximum magnetic quantum number. Expan- 
sion of Bj(fj,Bi/T) in fiHi/T at T > T N gives the linear 
temperature dependence for l/x(T,H) = H/M: 



X{T*) = T + T m (H) 
X (T,H) T* + T m (H) 

-a)M (J- 



(1) 



Here T m = (b-a)M (J+l)/6J, T» = 2T m , T m (B) = [1- 
cj(fiH/T) 2 ]T m , cj = (2J 2 + 2J + 1)/10J 2 , and x(T*) = 
9JT to /(J+1)mM , T m = 7T W where 7 = (6-o)/(6+o) < 
1. Eq. |T]) describes the data in Fig. 3 well, allowing us 
to estimate the density n of magnetic moments fi, which 
can provide the observed x- Taking xo = nfi 2 /3T m — 
1/395 from the inset of Fig. 3, T m = 12 K, and fx = 
3.62/ib for Nd 3+ 0, we obtain n = 3T mX o/^ 2 ^ 1.5 x 
10 22 cm -3 , close to the density of Nd atoms. Since the 
strong paramagnetism above T c is only consistent with 
an atomic density of magnetic moments fi ~ 3 — 4^b, 
Fe 3+ with much smaller fj, ~ 0.25 - 0.35/xs @, [3| can 
only add a few percent to x oc /i 2 . The crystalline field 



splits the degenerate orbital term of Nd 3+ 21, [22j into 
terms with spin moments /i smaller than [i of the free 
Nd 3+ but with n greater than the atomic density. 

We consider breakdown of AF by strong fields using 
a simple model, in which only spin degrees of freedom 
contribute to M, so the mean-field equations, 2M\ = 
Af tanh(^ffi/T) and 2M 2 = M tanh(^ff 2 /T) can be 
reduced to the following parametric equations 



z — cosh 1 [(2/8s) sinhs — coshs], 
m = smb. zj (cosh z + cosh s ) , 



I = 9s/2, h = z9/2 



(2) 
(3) 
(4) 



which define the dimensionless magnetization m = 
M/Mq and the order parameter / = L/Mq as functions 
of the dimensionless field h = /iH/Tj^ and temperature 
8 = T/Tff. Here the parameter s runs from to 2/9, and 
T N = M (b + a)/2. For T < T N , the AF order param- 
eter L(H,T) decreases as H increases, vanishing at the 
second order phase transition field H = H p : 

H p = [-fVT~e + e cosh- 1 (1 /Ve)]T N /n, 



(•5) 



at which M(H P ) = M \/l - 6. For H > H p , the magne- 
tization in the paramagnetic phase is described by 

m = tanh[(h-jx)/6] (6) 

Shown in the inset of Fig. 4 are m(h) and 1(h) curves 
calculated from Eqs. ©-© for high (6 = 5/11) and 
low (8 — 0.1) temperatures. At high temperatures m(h) 
has an upward curvature at H < H p and the downward 
curvature above the second order phase transition field 
Hp at which L(H) vanishes. At low temperatures, m(h) 
and 1(h) become multivalued, indicating hysteretic first 
order metamagnetic transitions due to spin flip in the 
sublattice with m antiparallel to H. 

One can immediately see a very different behavior 
of the mean-field M(H) as compared to the observed 
M P (H). Indeed, all M p (H) curves exhibit downward 
curvatures and can be described well by the Brillouin 
function M = M tarih((iH/T) with M = fin, n = 
(0.67 - 0.85 V B and n = (2 - 3) X 10 22 c m - 3 (the full set 
of fit parameters is given in the caption of Fig. 4). These 
values of n are about 2 times higher than the density of 
Nd 3+ in Nd(0,F)FeAs, which may indicate splitting of 
the 4 ig/2 ground term of Nd 3+ into different spin lev- 
els. Yet M P (H) shows no long-range AF behavior below 
T c , otherwise m(H) would be significantly depressed at 
H < H p , as depicted in the inset of Fig. 4. For exam- 
ple, if Tat = UK, 9 = 5/11, fi = O.SfiB, and 7 = 0.5, 
Eq. ([5]) yields the field H p ~ 16T required to destroy the 
AF state. However, M(H) in Fig. 4a exhibits neither 
the downward curvature nor metamagnetic transitions 
characteristic of the mean-field behavior at H < H p , 
while the data in Fig. 4b exhibit no sign of this behav- 
ior on a greater magnetic field scale either. Therefore, 
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the long-range AF ordering below T c is inconsistent with 
our data, in agreement with the conclusion drawn from 
neutron scattering measurements on La(0,F)AsF Q. At 
the same time, the Brillouin function for noninteracting 
spins gives a good description of the observed behavior 
of M(T,H) at T < T c , particularly the saturation of 
M (H) at higher H and decreasing the low-field slope of 
M (T) - n 2 nH/T as T increases. 

In polycrystals M(T, H) should be averaged with re- 
spect to easy magnetization directions in randomly ori- 
ented crystallites. For a uniaxial magnetic anisotropy 
and uniform distribution of easy magnetization axes in- 
clined by angles 6 with respect to H, we have, M = 
J* /2 M(iiH cos 0/T) cos0sin0d0 and X = X||/3 + 2xx/3, 
where X \\ = 2(1 - l 2 )T N /(a + b)[T + (1 - l 2 )~/T N ] and 
X±_ = \/b. If the in-plane AF structure is present, one 
can expect suppression of the AF state and local metam- 
agnetic transitions at lower T in the grains for which the 
in-plane field component H cos 9 exceeds H p . Also out-of- 
plane spin-flip transitions may occur if the perpendicular 
field component IT sin exceeds H c = \2K/(\j_ — Yii)] 1 ^ 2 , 



where K is the energy of magnetic anisotropy 2l|. 

Now we discuss possible mechanisms behind the ob- 
served effects. The first candidate would be paramag- 
netism of Nd 3+ ions whose exchange interaction is me- 
diated by conducting electrons on the AsFe planes. As 
follows from Fig. 3, this interaction could provide the 
AF ordering in the normal state below TV = TV/7 ~ 
10— 20K with 7 < 1. However, the superconducting tran- 
sition in the AsFe planes sandwiched between Nd(0,F) 
planes occurs at temperatures T c > TV effectively cut- 
ting off the 3D RKKY interaction (23| of Nd ions from 
different planes and reducing TV of the 2D XY moments 
of single Nd planes to zero. This may explain why the 
simplest Brillouin function describes our data so well. 
Since the paramagnetic signal is so strong that the den- 
sity of the magnetic moments n ~ (2 — 3) x 10 22 cm~ 3 is 
comparable to the atomic density of Nd ions, we can rule 
out that the behavior of M(T, H) is controlled by a small 
fraction of partially reacted magnetic second phases. 

Our high field measurements enable us to set a lower 
limit for TV of AF ordered Fe 3+ sublattice. In this case 
the longitudinal susceptibility x\\ at low T and H < H p 
is dominated by paramagnetism of Nd ions. However, as 
H becomes of the order of Tjq j \i the AF magnetization 
of the Fe sublattice becomes strongly nonlinear, exhibit- 
ing metamagnetic or spin flop transitions, which would 
manifest themselves in jumps on m{H) curves. Since no 
sign of such behavior has been observed up to H max = 33 
T, the long-range AF order can only be consistent with 
our results if H p given by Eq. (JSJ) is greater than H max , 
thus T/v > TV where H p (T p ) = H rnax . For H max = 33T, 
H = 0.35/zs |a 0, T = 20 K and 7 = 0.5, Eq. © gives 
T p « 21 K. However, much greater than 21K would 
be hard to reconcile with the Curie- Weiss behavior in 
Fig. 3 with T m = -/T N « 11 - 12 K, unless 7 < 1. 



Another possibility is a spin-glass state for which m(H) 
generally exhibits the behavior similar to that in the 
left inset in Fig. 4 with a broadened cusp around H p 
and possible kinks and hysteretic jumps [24]. Yet the 
behavior resembling that of M(H, 4.2K) has also been 
observed on some metallic spin glasses, although on a 
much smaller field scale (24j. Such a spin glass state 
could result from Nd ions and weakly coupled AF clus- 
ters of Fe 3+ , so that M = M | °° Ffai) tanh(^/T/T)d/i, 
where the distribution function F{^,i) can be expressed 
in terms of the distribution function of cluster sizes. This 
mechanism may also contribute to the fractional moment 
/i = J Q fiiF(ni)dfii obtained from the fit in Fig. 4. AF 
clusters of Fe 3+ around oxygen vacancies were reveal by 
the Mossbauer spectroscopy of SmOi-^F^FeAs @. 

In conclusion, doped Nd(0,F)FeAs can be regarded as 
a paramagnetic superconductor in which magnetization 
is dominated by the Curie- Weiss paramagnetism of the 
rare earth ions. Our high field magnetization measure- 
ments have shown no long-range antiferromagnetic order 
below T c despite noticeable AF correlations above T c . 
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